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We study the Fermi hquid properties of the cold atomic dipolar Fermi gases with the explicit dipo- 
lar anisotropy using perturbative approaches. Due to the explicit dipolar anisotropy, Fermi surfaces 
exhibit distortions of the d^2_-j^2-type in three dimensions and of the d^2_y2-type in two dimen- 
sions. The fermion self-energy, effective mass, and Fermi velocity develop the same anisotropy at the 
Hartree-Fock level proportional to the interaction strength. The Landau interaction parameters in 
the isotropic Fermi liquids become the tri-diagonal Landau interaction matrices in the dipolar Fermi 
liquids which renormalize thermodynamic susceptibilities. With large dipolar interaction strength, 
the Fermi surface collapses along directions perpendicular to the dipole orientation. The dynamic 
collective zero sound modes exhibit an anisotropic dispersion with the largest sound velocity prop- 
agating along the polar directions. Similarly, the longitudinal p-wave channel spin mode becomes a 
propagating mode with an anisotropic dispersion in multi-component dipolar systems. 



PACS numbers: 03.75.Ss, 05.30.Fk, 75.80.+q, 71.10.Ay 



I. INTRODUCTION 

The Fermi liquid theory is one of the most important 
paradigms in contemporary condensed matter physics, 
which sets up the fundamental framework to understand 
interacting fermion systemsii^. The most prominent fea- 
ture of the Fermi liquid theory is the existence of the 
Fermi surface and the long-lived low energy fermionic 
quasi-particles around the Fermi surface. The interaction 
effects, which may not be weak, can be conveniently de- 
scribed by a set of phenomenological Landau interaction 
parameters Fi in different partial wave channels char- 
acterized by the orbital angular momentum number I. 
Physical susceptibilities, including compressibility, spe- 
cific heat, and spin susceptibility receive important but 
still finite renormalizations by the Landau interactions. 
Furthermore, Fermi liquid states possess collective exci- 
tations such as the zero sound mode whose restoring force 
is provided by Landau interactions. 

When Landau interaction parameters are negative and 
large enough, Fermi liquid states may become unsta- 
ble toward Fermi surface distortions named the Pomer- 
anchuk instabilities'^ . The most familiar example is ferro- 
magnetism lying in the s-wave spin channel. High partial 
wave channel instabilities in both density and spin chan- 
nels have been studied extensively in recent years^"— . 
Another well-known, but of little practical relevance, 
is the Kohn-Luttinger superconductivity driven by in- 
teractions in the high angular momentum channels^^. 
The density Pomeranchuk instabilities in high partial 
channels exhibit anisotropic Fermi surface distortions 
which are electronic version of the nematic liquid crystal 
states^^— . Non-s-wave spin channel Pomeranchuk insta- 
bilities are essentially "unconventional magnetism" , in 
analogy to unconventional superconductivity. They in- 
clude both isotropic and anisotropic phases dubbed a 
and /3-phases as the counterparts of '^He-B (isotropic) 



and A (anisotropic) phases^—, respectivelj^. 

On the other hand, the rapid experimental progress 
of cold atomic physics provides an exciting opportunity 
to study quantum many-body systems with electric and 
magnetic dipolar interactions^"—. When electric dipole 
moments are aligned by the external field, dipolar inter- 
action decays as and thus is long-ranged in three 
dimensions (3D) but remains short-ranged in two dimen- 
sions (2D). More importantly, the most prominent fea- 
ture of dipolar interaction is its spatial anisotropy which 
possess the dz2-3r^ anisotropy in 3D and the dx^-y^ 
anisotropy in 2D, respectively. Considerable progress has 
been made in studying exotic properties in anisotropic 

condensations of the dipolar bosonic systems as reviewed 
inRef£i8-2i, 

Furthermore, dipolar fermionic systems provide an- 
other exciting opportunity to study exotic anisotropic 
many-body physics of fermions, including anisotropic 
Fermi liquid states and Cooper paring states. Physi- 
cal observable should exhibit the same anisotropy ac- 
cordingly such as the shape of the Fermi surface^^i'^ 
Anisotropic Cooper pairing and Wigner crystalliza- 
tion of dipolar Fermi gases has been theoretically 
investigated^^"— . Recently, Fregoso et aL™ studied the 
biaxial nematic instability in the dipolar Fermi gases as 
the d-wave channel Pomeranchuk instability, and general- 
ized the Landau interaction parameters to the tridiagonal 
Landau interaction matrices. 

Anticipating a great deal of experimental and the- 
oretical activity in the near future in polar molecular 
and atomic interacting fermionic systems, we have pro- 
vided in this article a comprehensive Landau Fermi liq- 
uid theory for polar interactions including the full effects 
of anisotropy. The explicit dimensionless perturbation 
parameter is the ratio between the characteristic dipo- 
lar interaction energy and the Fermi energy. The stan- 
dard textbook isotropic Fermi liquid theory is generalized 
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into the anisotropic version which exhibits many different 
features of both single-fermion and collective properties. 
Our theory is a leading-order perturbative theory in the 
polar interaction coupling constant, which is equivalent 
to a Hartree-Fock approximation of the interaction. We 
expect our theory to be quantitatively accurate in the 
weak coupling regime, but the qualitative aspects of our 
theory, e.g. the effect of anisotropy on the Fermi liquid 
parameters, should be generally valid. Our work only 
considers the explicit anisotropy due to the alignment of 
the molecular electric dipolar moment. The spontaneous 
anisotropic phase as a ferro-nematic state of the dilute 
magnetic dipolar Fermi gases has been recently studied 
by Fregoso et al. in Refj^ 

In Sect, nil the Fourier transforms of the dipolar in- 
teractions in 3D and 2D are summarized, which nicely 
exhibit the dr2_3z2 symmetry and the (1^2 _y2 symme- 
try, respectively. The anisotropic interactions result in 
anisotropic Fermi surfaces. The angular dependence of 
Fermi wavevectors, Fermi velocities, and effective masses 
are presented in Sect. IIIII All of them exhibit the 
same isotropy proportional to the the dipolar interaction 
strength at the leading order. 

Because the dipolar interaction mixes different par- 
tial wave channels, the Landau interaction parameters 
in isotropic systems becomes Landau interaction matri- 
ces. In Sect. IIVI the Landau interaction matrices in 3D 
calculated in Ref.'^S are reviewed, and those in 2D are 
constructed Fortunately due to the d22_3^2-symmetry of 
the dipolar interaction, the mixing turns out to result in 
the tri-diagonal matrix, between I and Z ± 2 in 3D. The 
thermodynamic susceptibilities would be renormalized by 
the generalized Landau matrices as will be discussed in 
Sect. lYl In this section, we will also explore the thermo- 
dynamic instability, where Fermi surface collapses along 
directions perpendicular to the dipole orientation. 

In Sect. IVIl we study the 3D collective modes in both 
the density and spin channels, respectively, focusing on 
the anisotropy effect. The dynamic collective zero sound 
mode exhibits an anisotropic dispersion. The sound ve- 
locity is largest if the propagation direction is along the 
north and south poles, and becomes softened as the prop- 
agation deviates from them. For directions close to the 
equator of the Fermi surface, the zero sound cannot prop- 
agate. The spin-channel collective modes do not exist in 
the s-wave spin wave channel. Nevertheless, well-defined 
propagating collective modes appears in the longitudi- 
nal p-wave spin channel, which has not been observed in 
condensed matter systems before. 

Conclusions and outlooks are made in Sect. IVIIl 



II. THE DIPOLAR INTERACTION IN THREE 
AND TWO DIMENSIONS 

The most prominent feature of the dipolar interaction 
is its spatial anisotropy, i.e., it depends on not only the 
distance between two dipoles but also the angles between 
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FIG. 1: Schematic sketch of a 2D dipolar system. The exter- 
nal electric field E controls the dipole orientation and thus 
the dipolar interaction. 

their relative vectors and dipole orientations. In the case 
that all the dipoles are orientated along the external elec- 
tric field E set as the z-axis, the dipolar interaction be- 
tween dipoles at ri and r2 reads 

(f 

V3D{ri-r2) = ^(1 - 3cos2%_fJ 

|ri - r2\-^ 

2cP 

= --p ^P2(C0S%_?J, (1) 

\ri - r2\-^ 

where Op^^p^ is the angle between ri — r2 and the £'-field; 
d is the electric dipole moment. The anisotropy exhibits 
in the angular dependence of the V^d with the form of 
the second Legendre polynomial. V^d is repulsive for 
9q < < TT — 00, and attractive otherwise, where 

00 = cos-^ ^ w 55°. (2) 
V3 

If the spatial locations of dipoles are confined in a two 
dimensional plane, and the £^-field is set in the a:z-plane 
with an angle of Oq relative to the plane as depicted in 
Fig. m then Eq. [T] reduces to 

(f 

V2D in - f2, 0o) = p ^ (1 - 3 sin2 9o cos^ 0) 

\ri ~ r2\'^ 

ff 3 

= "F ^{/2(cos6'o) - -sin^6'ocos2(/i}, (3) 

|ri - r2|'^ ' 2 

where is the azimuthal angle relative to the a;-axis. 
Eq. 13] is decomposed into the isotropic component and 
an anisotropic d-wave component, whose relative weight 
is tunable by varying the parameter angle ^o- If -E is 
perpendicular to the plane, i.e. Oq = 0, V2D is isotropic 
and repulsive. On the other hand, at Oq — Oq, V2D is 
purely anisotropic with the d-wave form factor cos 20. 
As the external electric field is tilted, i.e., 6q varies from 
to ^ , the 2D dipolar interaction gradually changes from 
an isotropic repulsive to an attractive one at large value 
of 9o. 

The anisotropics in the Fourier transform of dipolar 
interactions in both 3D and 2D exhibit in their depen- 
dencies on momentum orientations. In order to handle 
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both the long range part and the short distance diver- 
gence of the dipolar interactions, we introduce a small 
distance cutoff e beyond which the dipolar interactions 
Eq. [T]and Eq. [3] are valid, and a large distance cutoff R. 
e can be chosen, say, one order larger than the size of the 
dipolar molecule and R is the radius of the system. 



A. Fourier transform of the 3D dipolar interaction 



where Jo, 1,2(2:^) are the Bessel functions of the 1st, 2nd, 
and 3rd orders, respectively; I2 is the integral defined as 



hiqe^qR) 



qB. 



dxJ2{x), 



(6) 



and 72(0, +00) = 1. In the regular limit of qe — >■ and 
qR +00, the complicated form of Eq. [5] can be simpli- 
fied into 



For the 3D dipolar interaction, its Fourier transform 
can be performed as 

ji(9e) 



VsDiq) = Snd^ 



^^|F2(cos%) 



P2(COS0,-) 



(4) 



at ge — > and qR — 00. Its angular dependence is 
the second Legendre polynomial of momentum direction 
inherited from the real space form of Eq. [TJ ji (x) is the 
first order spherical Bessel function with the asymptotic 
behavior ji{x) — | as x — > 0; and ji{x) — ^ i sin(a:; — ^) 
as a; — >■ 00. 

Eq. m does not depend on the magnitude of q after 
the limits of — >■ and qR — !■ +00 are taken because of 
the spatial integral over renders the result dimen- 
sionless. At g = 0, Vsniq = 0) = arising from the 
fact that the spatial average of the 3D dipolar interac- 
tion is zero. However, V^Diq) is not analytic as q — >■ in 
the limit of i? — ^ -l-oo due to the angular variation. In 
fact, the smallest value of qR is at the order of 0(1), this 
non-analyticity even exist for large but finite value of R. 

An interesting feature of the above Fourier transforms 
is that Vsoiq) is most positive when g is along the z-axis, 
but most negative when q is in the equator plane, which 
is just opposite to the case of Vsoir) in real space. This 
can be intuitively understood as follows. Considering 
a dipole density wave propagating along the z-axis, the 
wave fronts (equal phase lines) are perpendicular to the 
dipole orientation, thus the interaction energy is repul- 
sive. On the other hand, if the dipole density wavevector 
lies in the equator, then the wave fronts are parallel to 
the dipole orientation which renders the interaction neg- 
ative. 



B. Fourier transform of the 2D dipolar interaction 

In two dimensions, the Fourier transform of Eq. |3] is 
more subtle, which can be expressed as 

I e R 



q[Jl{qe) - MqR)] - ql2{qe,qR)} 



ird sin 9, 



'ocos20g gj; 



-J2{qR) J2{qe)- 



qR qe 
[UqR)-Mqe)]+l2iq^,qR)}, (5) 



V2D{q,Oo) = 2Trd'P2{cos9o){--q) 

+ nd^ sin^ 9o qcos2(j)q. (7) 
In particular, at Bq — 9q, Eq. [7] is entirely anisotropic. 



2'Kd^ 

V2D{q,0o = 9o) = ——qcos2(j)q, 



(8) 



as one can see from the real space interaction in Eq. |3l 
Eq. [7]even holds in the long wavelength limit qR 0(1). 
On the other hand, for the short wavelength limit x — qe 
at the order of 0(1), additional terms should be added 
into Eq. [7] as 

AV2D{q~^0,9o) = 2Trd^P2{cos9o)-[xMx) + xl2{0,x)] 

e 

+ nd'^ sin^ 9q cos (f)q-[3J2{x) + xJ:i{x) 



xhiO, x)] 



(9) 



In the dilute limit, where the Fermi wavevector kf satis- 
fies kfe <^ 1, we only need to use q up to the order of fcj. 
Thus in the sections below, we shall use Eq. [7] as the 2D 
Fourier transform of our dipolar interaction. 

The anisotropic component of V2D{q) (second term of 
Eq. [7]) has a similar feature to the 3D case, which is 
repulsive along the x-axis but negative along y-axis. No- 
tice that there is also an isotropic component for V2D 
which can be either positive or negative depending on the 
external parameter 9o. Therefore, in momentum space, 
the dipolar interaction in 2D contains both s and d-wave 
components, while it only has a d-wave symmetry in 3D. 



III. HARTREE-FOCK SELF-ENERGY AND 
FERMI SURFACE DEFORMATION 



The anisotropy of the dipolar interaction exhibits in 
the fermion self-energy at the Hartree-Fock level, which 
results in anisotropic Fermi surface distortions as studied 
in Ref o^^i^^'^° . A convenient dimensionless parameter to 
describe the interaction strength is defined in 3D and 2D 
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pj3D 
kfo 



d'^mk^P 



2D 



^int Jo 



K/O 



(10) 
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FIG. 2: The dipole moments polarized by the external E- 
field are along the z-axis. The polar and azimuthal angles of 
k and k' are denoted by ^^(fc') and 4>k(k'), respectively. 0^_j:, 
is the polar angle of the vector difference oik — k' but not the 
angle between k and k' . 
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FIG. 3: The dimensionless fc-dependent function Izd(y^) of 
the 3D Hartree-Fock self-energy in Eq. 1131 



A. Anisotropic Hartree-Fock self-energy in 3D 

We first consider the situation in 3D, the Hartree-Fock 
self-energy can be expressed as 

= ^^Z^^2(cOS^?^_£,)^fc' 
k' 

where is the Fermi occupation number; the Hartree 
contribution goes to zero. In Eq. [TTl the leading con- 
tribution of the anisotropy to I]^^(fc) comes from the 
interaction. The effect from the Fermi surface distortion 
is at a higher order of A, and thus is neglected, thus we 
will take the Fermi wavevector in Eq. [TT] as that of zero 
dipolar interaction kj^. dj:_g is the angle between the 

momentum transfer k — k' and E, which satisfies 

(k-k')-E 

\k-k'\ ^ ' 

Please notice that is not the angle between k and 

k' (see Fig. 

This Hartree-Fock self-energy can be evaluated analyt- 
ically as 

Ef/(fc) - ~2\Elf^P2{cos6u) I3d{^), (13) 

fo 

where Ek^g = h'^k^f^/2m and 



where fc^g ' are the non-interacting Fermi wavevectors 
in 3D and 2D, respectively; E^t is the average dipo- 

lar interaction and Ei, ' — tt are the kinetic 

energies at the Fermi surface in the absence of interac- 
tion. Since that X^d and X2D are at the same order, we 
will not distinguish them but use A = X^d for both 2D 
and 3D. In this section, we shall evaluate the Hartree- 
Fock self-energy in 2D and 3D perturbatively at the lin- 
ear order of A from which the Fermi surface distortion 
will be determined. We will see that both the Hartree- 
Fock self-energies and Fermi surface distortions possess 
the same symmetries inherited from the corresponding 
interactions. 



vr f , 3 3(l-a;2)3 ,l + a;,i 

(14) 

is a monotonically increasing function as depicted in Fig. 
[3l This form agrees with the self-energy calculated for 
the 3D dipolar Fermi gases in Refi^. Right on the Fermi 
surface, 

A-TT 

Sf/(fc = fc??) = --XElf^P^icosek). (15) 

Naturally Eq. [13] exhibits the d-wave form factor 
P2(cos 9k)- We shall see more examples below where the 
physical quantities possess the symmetry originated from 
the dipolar nature. 
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2.0 k. 



FIG. 4: The fc-dependent components of the 2D HF self- 
energies in 2D. The sohd (blue) and dashed (red) curves rep- 
resent the dimensionless functions lio'ik/kf) and l2D{k/kf) 
defined in Eq. [16] 



B. Anisotropic Hartree-Fock self-energy in 2D 

Similarly, for 2D system, the HF self-energy is evalu- 
ated as 



^2D 



h,2D 
'«/0 



- sin^ 00 cos 2(j,k /2T (t2u) I ' (16) 

where fc^^ is the Fermi wavevector at zero dipolar inter- 
action and 

lif^{x) = 127r^ dx'x'\x + x'\E(^^^'^ 

The functions K{y) and E{y) are the standard complete 
elliptic integral of the first and second kinds, respectively 

as 



K{v) 

E{y) 



da- 



1 



— j/^ sin^ a 
day 1 — j/^ sin^ a. 



(18) 



The behavior of /2d(^) ^-nd I2d^{x) are plotted in Fig. 
m Eq. [TBI indicates that the self-energy in 2D comprises 
of an isotropic and an anisotropic terms. The former 
shifts the chemical potential ^, while the latter distorts 
the Fermi surface. At fc = kjQ, Eq. [16] reduces into 

Sf/(fc = fcf?,^fc,^o) 
= -XElf^^i-^ P^icosOo) + i^sin2 0ocos20fe}. 

(19) 



Again, Eq. [16] contains the same angular dependence as 
the 2D interaction does. Both Eq. [T3] and [TH] indicate 
that the Hartree-Fock-self energy changes monotonically 
with the momentum k in 2D and 3D (Fig. [3] and Fig. 



C. Hartree-Fock self-energy for two-component 
dipolar systems 

In two-component dipolar systems, the exchange 
part of the Hartree-Fock self-energy only exists intra- 
component interaction, and the Hartree part exists in 
both intra and inter-component interactions. Since the 
Hartree contribution vanishes in the 3D dipolar systems, 
the Hartree-Fock self-energy remains the same as in the 
single-component systems. In the 2D dipolar systems, 
there is an extra contribution from the inter-component 
Hartree interaction as 



AyHF _ 



3tt^ 

--^XEkfoP2{cos 9o), 



(20) 



where e is the 2D phenomenological cutoff constant we 
discussed before. This extra self-energy contribution is 
momentum independent and can be offset by an overall 
shift of the chemical potential. 



D. Anisotropic Fermi surface distortions 

The anisotropic Hartree-Fock self-energy from the 
dipolar interaction naturally results in anisotropic Fermi 
surface distortions. With Eq. [T3] and [TB] we determine 
the distortions by solving: 



eHF{kf) = eo(fc/) 



(kf) = ^iin,X). 



(21) 



Here, n is the particle density and we recall that kj^ is 
the Fermi wavevector for A = 0. With nonzero A, the 
dependence oi kf on the polar angle in 3D is solved as 

kfiOk) 



l,3D 
'«/0 



l-—X^ + —XP2{cos0k). 



(22) 



Similarly, in 2D, the dependence of k on the azimuthal 
angle is solved as 



h.2D 



l_l??!siu^0oA^ 



25 



sin 6'oAcos2(/)fc 



(23) 



The anisotropic Fermi surface distortions in 2D and 3D 
show linear dependence on dipolar interaction strength 
A. Due to particle number conservation, there are small 
shrinking of Fermi momentum in both 2D and 3D at the 
quadratic order of A. Notice that these two equations 
are correct to the order of A. The A^ terms appear to 




FIG. 6: The 2D distorted Fermi surfaces for A = 0.02. The 
soUd (black), dashed (red) and blue (dotted) curves corre- 
sponds to an external electric field configuration at 6^0=0, Oq 
and |- respectively. The electric field is projected to x-axis in 



FIG. 5: The deformed Fermi surface for X = ^ in the 3D ^^^^ system, 
dipolar system obtained by the perturbative (solid) and vari- 
ational (dashed red) approaches. The external electric field 
lies along the z a^cis here. 



conserve particle numbers due to this lowest order Fermi 
surface deformation. When higher order contribution is 
considered, anisotropy can also enter in the second order 
of A and particle conservation will in turn lead to an 
isotropic correction. 

We emphasize here again that the results obtained 
based on Eq. [13] and Eq. [16] are perturbative and are 
correct for A <C 1, while they provide qualitative features 
for A < 1. As a comparison with the variational ap- 
proach suggested in Refi^, we plot the deformed Fermi 
surface in Fig. [S]for A = ^ in 3D, which corresponds 
physically to the stability limit of the dipolar system (to 
be discussed in Sect. [V]). Our perturbation results give 
rise a less prolate shape of the Fermi surface than those 
obtained from the variational approach. Fig[B] shows the 
same distorted Fermi surfaces in 2D. At 6*0 = 0, the in- 
teraction is isotropic and does not deform the Fermi sur- 
face. It becomes prolate when 9q increases. We would 
like to emphasize that the interaction is most anisotropic 
at 6*0 = ^Oi though a larger 6*0 will stretch the Fermi 
surface more. 



E. Anisotropic Fermi velocities and effective 
masses 

The Fermi velocities become anisotropic in both 2D 
and 3D at non-zero A. In comparison to interacting elec- 
tron systems, the Hartree-Fock contributions to the 2D 
and 3D dipolar Fermi gases are less singular, i.e.,, the 
radial derivatives at the Fermi surface remain continu- 
ous: the Hartree correction of Fermi velocity is continu- 
ous because the Fourier transforms of the dipolar inter- 
action are finite as g — > 0. Furthermore, the difference of 
the Hartree-Fock self-energies for the single component 



and multi-component dipolar gases is momentum inde- 
pendent, thus the Fermi velocities are the same for both 
cases and we will not distinguish them in this subsection. 

We first consider the 3D case. After taking into ac- 
count the Hartree-Fock correction, the projections of 
Fermi velocities along the radial direction k and angu- 



lar direction 

as 



are defined as k 



v^^{k) and eg 



-3D 



(k) 



773Dt 



f)yHF 



dk 



ayHF 



{k,ek 



(24) 



Taking into account the angular dependence of Fermi 
wavevector in Eq. 1221 and we arrive at the radial and 
angular Fermi velocities corrections to the linear order of 
A as 



l,3D 



1.3D 



{l- |AP2(cOS0fe)}, 



ttA sin 29 k. 



(25) 



(1) the 



There are two opposite contributions to fc/ 
Fermi surface deformation, kf{9k) and (2) the Hartree- 
Fock self-energy modification. In 3D dipolar system, the 
latter is stronger and thus the radial Fermi velocity is 
suppressed at north and south poles but enhanced along 
the equator of the Fermi surface. 

According to Eq. [25] we define the longitudinal and 
transverse effective masses m 



3D, 



and m'^jj j_ in 3D as 



1 



1 



kj^'iOk) 
k^f^iOk) 



-|l-7rAP2(cos0fc)|, 

771 L J 



= — 7rAsin26'fc. (26) 

777, 
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For the 2D case, the projections of Fermi velocity in 2D 
along the radial k and the azimuthal angle e^j. -directions 
are calculated as 



U2D 



,1 + A 

T 



Asin^6'osin2(/)fe. (28) 



For a general value of 7^ ^o, there is an additional 
isotropic renormalization to v'j^ {cfyk, 9q) in Eq. 1271 which 
comes from the isotropic part of the 2D dipolar interac- 
tion. 

This means that in 2D, we also have a reduced radial 
Fermi velocity at (pk — 0, tt, while it is boosted when (jj^ = 
Accordingly, similar to the 3D situation, we can also 
compute the longitudinal and transverse effective masses 
rri^j-) II and rn^j^^x^ in 2D as 



interaction at the Hartree-Fock level. We define the 3D 
differential DOS of the single component N^^{Vlk) in the 
direction of ilk as 



|H- A[47rP2(cos6lo) 
sin2 0ocos20fe]}, (27) N^'^i^k) 



k}D{dk)dnk 

(27r)3w3^( 

j-^[l + —XP^{cos9k)]dnk. (30) 



At the linear order of A, N^-^{flk) develops the same 
anisotropy of P2(cos6'). The total DOS reads 



3D 



dflk 

An 



(31) 



= l|l + A 
m L 

1 



47rP2 (cos ^^o) sin'^ cos 2(/)fc 

5 



1 IGtt 



'^2D ±i4'k,9o) m 5 



A sin 9osin2(f>k- 



(29) 



F. Renormalization of density of states 

Now we study the renormalization of the density of 
states (DOS) at the Fermi surface due to the dipolar 



which does not change at the linear order of A compared 
to that of the free Fermi gas. This means that the specific 
heat, which is proportional to the total DOS at the Fermi 
surface, is not renormalized to the linear order of A. In 
considering the actual area of Fermi surface is enlarged, 
the correction at higher orders of A should increase the 
total DOS. 

Similarly, we consider the anisotropic case in 2D. The 
2D differential DOS of the single component N^^{^) at 
the Fermi surface along the direction of the azimuthal 
angle 0^ is defined as 



27r 



[fcr(<^fc'^o)]2 + ( 



de.D 
-J—Y 

d(j>k 



m,d(f>^ 



d(t)k f 
7r)2 I 



1 + A 



47rP2(cose'o) 



147r 



■ sin 9o cos 20^ 



(32) 



which clearly exhibits the d-wave anisotropy. The in- 
tegrated total DOS at 2D is not changed at the linear 
order of A either. Similarly, the 2D specific heat is not 
renormalized to the linear order of A. 



IV. LANDAU INTERACTIONS FOR THE 
DIPOLAR FERMI LIQUIDS 

In this section, we construct the Landau Fermi liquid 
Hamiltonian for the dipolar fermion gas which has also 
been studied by Fregoso et. before. In the isotropic 
systems, the interaction effects in the Fermi liquid the- 



ory are captured by a set of Landau parameters Fi in 
different partial wave channels. In dipolar systems, the 
anisotropic interaction leads to the mixing of interactions 
in different partial wave channels, thus we need to gener- 
alize the concept of Landau parameters into the Landau 
matrices. 



A. Landau interaction matrix for the single 
component dipolar gases 

In order to have a common reference Fermi surface for 
different but small values of A, we choose it as that of the 
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free fermion gas with the same particle density. Define 
the variation of the Fermi distribution at momentum k, 

Srij^ ^ - no^k, (33) 

where no(fc) = 1 — 0(fc— fc^^). Assuming a small variation 
of the fermion distribution Snj: close to the Fermi surface, 
the ground state energy of the Fermi liquid state changes 
as 

SE^J2 ^fc^^fe + ^ E ^(^' k')Sn^Sn^„ (34) 

k,k' 

where k, k' are momenta close to the Fermi surface; 
/(fc, k') is the interaction function describing the forward 
scattering amplitude, efc includes the bare parabolic dis- 
persion e° and the self-energy correction T,HFik) which 
exhibits the d^^-sz^ anisotropy. At the Hartree-Fock 
level, f{k, k') is expressed as /(fc, k') = V[q = 0) — V[k — 
k') where the first and second terms are the Hartree and 
Fock contributions, respectively. Due to the explicit rota- 
tional symmetry breaking, f(k,k') depends on directions 
of both k and A;', not just the relative angle between k 
and k' as in isotropic Fermi liquids. For the 3D dipo- 
lar system, the classic Hartree term vanishes because the 
spatial average of the dipolar interaction is zero. How- 
ever, due to the non-analyticity of the Vsoio) as (f — > 0, in 
the calculation of the zero sound collective excitation in 
Sect. IVI[ the dependence of q in the Landau interaction 
needs to be included as 

k,k' 

where /(fc, k'; q) = V{q) — V{k — k') and n-^ ^ = ct^^cj. 
1. 3D Landau interactions 



dipolar interaction, its spherical harmonics decomposi- 
tion becomes 



(38) 

where remains diagonal for m but couples partial 

wave channels with /' = 1,1 ±2. This is a direct result 
from Wigner-Eckart theorem because the dipolar inter- 
action possesses the symmetry of d^^^^z^. The mixing 
between and I with V = I ± \ la forbidden because the 
dipolar interaction is parity even. 

For the 3D systems, the matrix elements fw-^rn are 
tridiagonal and have been calculated by Fregoso et al. 
in Ref.^° as 

f?i^m - d'{a[llSu'+a\^Sij,^2 + a^lS,,.2) (39) 
where 

(1) _ A7t{P + I - 3m^){2l + 1) 
~ /(/-|-l)(2/-f 3)(2Z- 1) ' 

(2) ^ 27r 

{I + + 2){2l + 3) 

X ^[{l + l)^ ~m^][{l + 2f-m% (40) 

except for the I = I' = m = channel, where we have 
foo'oi'l) — ^3D{q)- Please note that in Eq. [351 we use 
the standard normalization convention in Ref.— which 
is different from that in Ref.^°, thus the parameters in Eq. 
1401 are modified accordingly. Sign errors in the original 
expressions of Refi^ are corrected here. It can be proved 
that for each I, /j^P^'s satisfy the relation that 



In this part, we review the Landau parameter calcula- 
tion performed in Ref."^°. Since quasi-particle excitations 
are close to the Fermi surface, we integrate out radial 
direction and obtain the angular distribution. In the 3D 

system, 

/k'^dk 
J^Sn^k- (36) 

We expand Sn{Q^) in terms of the spherical harmonics 
as 

Sn{nj:) = ^ Yijn{nj:)Sni^, (37) 

Im 

where l^m's satisfy the normalization convention 
/ dil Yi*^{fl)Yim{^) = 1. Due to the anisotropy of the 



E = 0- (41) 



We define the average effective radial mass as 

rn^D ^ ^ J c^^^fe'^sc.lK^fe)- (42) 

The perturbative result in Eq. [521 correct to the linear or- 
der of A shows fh* = m. Following the standard method 
to make the Landau matrix dimensionless, we multiply 
the single component DOS. Fjj^^ = SLn-^Djsd^ ^^lere 
N^D ^ (rnkj^)/{2hn^) is the DOS of free Fermi gas. 
To gain some intuition, we present some values of the 
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low order Landau matrix elements in terms of A as 



order as 



Fi^i^oi^,-) = 47rAP2(cos0,), F^^^q = -ttA; 
187r , _ 97r , 



r3D 



-^22:0 — ~J~ ' ^22;±1 " ^22;±2 — 

p3D _ 37r 3^ _ V&TT 

p3D _ l^TT _ A, p3D _ n 

-^^33:0 — ^^'^J-f^33:±l — ' 33;±2 — 'Jj 



p3D 



33;±3 



A, 



(43) 



^22 ^ 
-,2D 



127rP2(cos6'o)A; Fi-^i^ = -47rP2(cos 6lo)A; 
An 

-—P2{cos0o)X- 



Fq2 = -7rsin^6loA; F^^^i = Stt sin^ 6I0A 



(48) 



Eq. [Sni and Sni generalize the usual Landau parameters 
in isotropic Fermi liquid systems to Landau matrices in 
the dipolar systems. In fact, matrix formalism is a natu- 
ral generalization as long as anisotropy enters the system. 
It is the dipolar nature that constrains our Landau ma- 
trices to be tridiagonal as shown above. 



2. 2D Landau interactions 



B. Landau interaction matrix for two-component 
dipolar Fermi gases 



Similarly, in the 2D system, we define the angular dis- 
tribution 5n(0g) and its decomposition in the basis of 
azimuthal harmonics e™"^*^ as 

^n,, = ^e™*.-Jn,„. (44) 

The Landau interaction can be represented by a matrix 
as 



/2^(fc,fc') = fnZ'^^'""^'^' 



(45) 



where /^^/ is non-zero when m' ~ m, m ± 2. 

We further present our calculation for fmm' in 2D 
which reads 

/™™'(^o) = kj^d'{P2icoseo)bi^^S,„^^. 

+ sin2 0o(6(^)5,n,™'-2 +6|^^5rn',™-2)}, (46) 



where 



6(1) = - 



(2m- l)(2m-h 1)' 
2 

'(2m-f l)(2m-|-3)' 



(47) 



At 9o — do where P2(cos^o) — 0, the dipolar interaction 
is purely anisotropic (see Eq. [5]) and the diagonal matrix 
elements vanish. At this particular angle, the interaction 
has angular momentum m = 2, and the matrix element 
is non-zero only when m' = to ± 2 as a result of the 
Wigner-Eckart theorem. 

Similarly, for 2D, we also define the average radial ef- 
fective mass m* which equals to to at the first order of 
A. After multiplying the 2D DOS of the single compo- 
nent Fermi gas, the Landau matrix becomes dimension- 
less = ^^o'''/™™' where N^^ = m/{2TTh^) is the 

DOS of 2D free Fermi gas. Some low order Landau ma- 
trix elements are presented in terms of A at the linear 



We next consider the Landau interaction matrix a 
two-component dipolar Fermi gas in 3D. The intra- 
component Landau interaction contains both Hartree 
and Fock contributions as in Eq. [351 The inter- 
component Landau interaction only contains Hartree 
contribution. The general Landau interaction function is 
decomposed into density channel response fs and spin- 
channel response as fa 

fa^^sik, k') = r{k, k')6^p5^s + /"(fc, k'Wp(T^s, (49) 

where fs and fa at the Hartree-Fock level are expressed 
as 

r(fc,fc';g) = V{q)~]^V{k~k')- 



r{k,k';q) = --V{k-k'). 



(50) 



Because the DOS at the Fermi surface for the two- 
component Fermi gases is doubled compared to that of 
the single-component gases, the Landau matrix elements 
in the density channel F^f!^ and in the spin channel 

Ffi^m Hartree-Fock level equal to those F^R^^ de- 

fined for the single component case 



T^3D,s 



F, 



3D, a 



3D 



II' :m 7 



(51) 



if at least one of / and I' are nonzero. The case of I 
Z' = is special, for 3D we have 



FoOiO - 2PoO;0; ^00;0 " 0' 



(52) 



In 2D, a similar conclusion applies at the Hartree-Fock 
level as 



p2D,s _ p2D,a _ p2D 
mm' mm' ^ mm' ^ 



(53) 



if at least one of to and to' are nonzero. For to = to' = 0, 
we have 



^00""'^'" - F^n ± A^^P2(cos^?o), (54) 
ekj ' 



where e is the short range cutoff defined in Sect. |TT1 
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C. Landau interaction matrix for A'^-component 
dipolar Fermi gases 

The general A^-component case is essentially similar 
in which A'^ arises from the hyperfine multiplets. The 
SU{N) symmetry is very accurate since the electronic 
dipolar interaction is independent of the internal hy- 
perfine components. A Fermi liquid theory for the 4- 
component Fermi gas with S'J7(4) and S'p(4) symmetry 
has been constructed by one of us in Ref."^^''^'^, which can 
be easily generalized to the A'^-component case here. For 
the convenience of presentation, we first define our con- 
vention of the A*"^ — 1 generators of the SU{N) group 

[Ml-\k = SiiSjk + SikSji {l<i<j< n), 



\lk 



diag(l,...,l,-(i-l),0,...) 

Vj(j-1)/2 



(2<J <n), 



(55) 



where M\ and are the SU{N) version of the 
Pauh matrices of ci, a2 and 0-3, respectively. Then the 
SU{N) Fermi liquid Landau interaction can be written 
as 

+ ^m3„^m3^,}, (56) 

i 

where fs and fa are expressed as 

nk,k';q) = V{q)-^V{k-k'); 

nk,k';q) = -^Vik~k'). (57) 

Again due to the DOS at the Fermi surface is A^- 
times enhanced, the Landau matrix elements for the A^- 
component dipolar gas in the density channel F^!^^ and 
in the spin channel F^!^^ at the Hartree-Fock level equal 
to those FjjPjj^ defined for the single component case as 



IT'S Tpa 77 

^U';m — ^ll'-m — ^ II' \m 



(58) 



when at least one of I and /' are nonzero. When I — I' — 0, 
we have 



-^00:0 = NFoq-q; F^g g - 0. 



(59) 



In other words, Fqq.q has a large iV-enhancement com- 
pared to all of the other Landau matrix elements. 

Again, in 2D, a similar conclusion applies at the 
Hartree-Fock level as 



p2D,s _ p2D,a _ p2D 
mm' mm' mm' ' 



(60) 



if at least one of m and m' are nonzero. For m = m' = 0, 
we have 



F 



2D,. 
00 



^^o'o^-f A(Ar-l)-^P2(cos0o), 
ek^ ■ 

37r2 



^0 ■ - ^00 -A-2^P2(cosf?o). 
ekj ^ 



(61) 



V. THERMODYNAMIC QUANTITIES 

We next consider the renormalization to the thermo- 
dynamic properties from Landau interaction matrices. 
For the simplicity of presentation, we only consider the 
single-component dipolar systems here. With slight mod- 
ifications, the results also apply to the multi-component 
cases. 



A. Effective mass and Landau interaction matrix 
elements 

In this subsection, we rederive the effective mass renor- 
malization discussed in Sect. IIII El using Landau matrix 
formalism. We will see how the Landau parameters in 
Eq. [39] and |46] enter in the effective masses in 3D and 
2D respectively. The formalism below is general for any 
anisotropic Fermi liquid system (with azimuthal symme- 
try in 3D). 

In Galilean invariant systems, the fermion effective 
mass renormalization ra* /m = (1 -|- ^i^f*) is an important 
result of the isotropic Fermi liquid theory. This results 
in the same renormalization factor for the DOS at Fermi 
surface and the specific heat as Cfl/Cfg = 'm*/m, 
where Cfl and Cfg £^re specific heat for Fermi liquid 
and ideal Fermi gas, respectively. 

For the anisotropic 3D dipolar systems with the 
Galilean invariance, we present a similar result. The re- 
lation connecting effective mass and bare mass still holds 
for anisotropic interactions2ii2^ 



deik) 
dk 



k_ 

m 



rf'fc'.^3D(^.^,^5n(6(fc')) 



(27r)3 



dk' 



(62) 



However, for anisotropic systems, a self-consistent solu- 
tion to Eq[62]has to be done numerically. To the linear 
order of A, we perform the analytic calculation by ap- 
proximating e(fc') in the RHS of Eq. [55] with the free 
fermion energy as follows. We take the radial derivative 
of Eq. [62] 



1 



m 



3D, 



1 

m 
1 

m 



dnk'f"{k-k'){k-k') 



1 - Fii^\\{Ok) - Fi3 \\{9k) 



(63) 



where 

m F, 



3D, 



9k) is the radial effective mass given in Eq. 
and -Fi3,||(6'fc) in Eq. [33] are the angular 
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dependent Landau parameters defined as follows 

= ^P2(cos0fc), (64) 



We define the variation of the angular distribution is de- 
fined as 



f P 

5n{Qk) = j dk-^^^—^5n{k,Vlk), 



(71) 



-^13,11 (^fc) 



4 F^^ 

- J2 ^Y*^{e,,o)YUOk,o) 



m=0,±l 

2T7r , 



3 , V21 

-AF2(cos dk). 



where only the radial integral of k is performed. The total 
density variation can be expressed a.s Sn = J dftkSn{nk)- 
The spherical harmonic expansion of (5n(f2fc) is defined 
as 



(65) 



(72) 



Eq. IMlis a generalization to that of the isotropic case of ^j^^ corresponding ground state energy is represented 



m* /m — 1 + which can be rewritten as 

1 ^ 1 iVg^/f 
m* m m 3 

Thus to the linear order of A, 

^ ^ {l-^AP2(cOS0fe)}, 



as 



m*J6k) m 



(66) 



(67) 



(73) 



which agrees with Eq. [511 At the linear order of A, 
as an averaged effect, there should be no specific heat 
renormalization due to the dipolar interaction. 

Parallel analysis can be carried out for the 3D trans- 
verse effective mass as 



where the first two terms are the variation of kinetic and 
interaction energies respectively, and the last term is the 
coupling to the external fields with partial wave channels 
of Im. Expanding the Hartree-Fock single particle energy 
around fcj^ as 



e^j^{k,n,) = eo{k}^) + j:^j^{kjf^,n,) 



3D\ I s^HFn„3D 



3D 



1 



1 



rf=i //) \ , 7=1 //I M ^ro\ The variation of the kinetic energy is represented as 



"''3D,± 

where 



V 



= / dl7,^2^[ '"^^'ll^"'V o^^]-Mfa(»fc)]^ 



gi^ll;0 
47r f2 



) sin 29k 



3 V 21 



Fli'^^Y^i{eu,o)Y3a{ek,G) 



(y3,~i(^fc,o) 



F, 



•3D 



13;-1 



^ ^^^\k)^,nk)8n{nk)], (75) 

Eq. [75] can be expressed as 
5Ek, 



V 



(69) 



(76) 



For the 2D effective mass, the Landau parameters renor- 
malizations are summarized below: 

L 0. - 1[1-F--(F,- +F-)cos2</>..], 



where /i^q = |i/|Ai^|^, and to* is defined in Eq. |32] 
The perturbative results at the linear order of A give rise 
to 



M'^^m = ^'illSw + m\:^Su'-2 + m\:>Su,^2, (77) 



"^^20 ±i^k,Oo) TO 



where 



(70) 



They all agree with the previous results in Sect. IIII El 
based on Hartree-Fock calculation, just as expected. 



B. Thermodynamic susceptibilities 

We study the variation of the ground state energy of 
the dipolar Fermi gas respect to Fermi surface distortion. 



„ _ (P + 2^-3to^) 

'™ " ^ (2Z + 3)(2/-l) 
_ + (1) 

" ^4(2^ + 1) ' 



(78) 



(2) 



3tt 



'[(/+ 1)2 -to2][(; + 2)2 -to2] 



2(2? + 3) y (2;-f l)(2Z-f 5) 

3{l + l){l + 2) (2) 

u> „ A. 



Ay'{2l + l){2l + 5) 



(79) 
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The variation of the interaction energy reads 

SEint 1 rSD r * jr 

= 9 2^ Jw ,,Jnimdnim- 



V 



(80) 



Im 



With the total field him = /ifm + ^20 ' ^^'^ ground state 
energy is represented as 



SE 



= SEkin + 5Eint - 4:TThimSnir. 



= 47r 



/r ;m 
m* at3D 



(81) 



where xo — ■ The matrix kernel Kw reads 



3Z? 



^(2; + l)(2Z' + l) 



(82) 



The the expectation value of the 6ni„i in the field of him 
can be straightforwardly calculated as 



5nim = Xi){Kfi^,m) 



(83) 



Thus XoiKfi^m)^^ ^^'^ renormalized susceptibility ma- 
trix for a 3D dipolar Fermi system. 

A parallel study can be applied to 2D, where we re- 
place the quasiparticle density fluctuation Sni„i by Sum 
as defined in Eq. 1441 The corresponding result is similar 
to the 3D case, except we replace the M^^^ matrix by 



K. 



2D 



M. 



2D 



F 



2D 



(84) 



where 



M. 



2D 



1 + 47rAP2(cos6'o) 



— A Sin^ 6*0 (Sm.rn'-2 + Sm' .m-2) ■ (85) 

5 

The renormalized susceptibility in 2D has the same 
form as in Eq. [551 



C. Thermodynamic stability 

In isotropic Fermi liquids, Fermi surface becomes un- 
stable if anyone of the Landau interaction parameters Fi 
is negatively large enough, i.e., F; < —{21 + 1). This can 
be understood by treating Fermi surface as elastic mem- 
brane. The kinetic energy always contributes to positive 
surface tension, while interaction contributions can be 
either positive or negative depending of the sign of Fi 
in each channel. If the negative contribution from in- 
teraction exceeds the kinetic energy cost, Fermi surface 
distortion occurs. 

In the 3D anisotropic dipolar Fermi gas, we diagonalize 
the interaction matrix Kfi^^ as 



r^3D 

J^ll'ry 



Tm'diag{fl^\t,T,t^"2\-}Trr 




N 0.4 



FIG. 7: A plot of the lowest eigenvalues /i^' of Kfi?.^ versus A 
for m = 0, ±2. The instabilities of the s-channel (solid blue) 
and d±2-channel (dashed red) are revealed by the vanishing 
of their eigenvalues when A passes 0.135 and 0.35 respectively. 



The thermodynamic stability conditions can be similarly 
stated as each of /x™ is positive, i.e., 



> 



(87) 



for arbitrary m and i. It is not difficult to observe that 
in 2D, we simply need to replace the matrix by K'^^, 
defined in Eq. [M] In the isotropic systems, KfiP.^ be- 
comes diagonal, and this stability criterion reduces back 
to that of the Pomeranchuk. 

For our anisotropic dipolar system, we diagonalize the 
if-matrix numerically and determine the instability con- 
ditions. For example, we compare two strongest insta- 
bilities in the sectors of m = and that of to = ±2 
with even I. The former mainly lies in the s-channel and 
the latter mainly lies in the rf3.2_j^2^2i2;y"Channel, both of 
which hybridize with other even partial wave channels 
with the same values of m = 0, ±2, respectively. Their 
eigenvalues are denoted as /Zs and ^id±2 1 respectively. The 
i^Q(^g((?) explicitly depends on the orientation of q. We 
put q along the equator which results in the minimal 
eigenvalues of Kfi^, and plot /ig and fid±2 Fig. [7| 
The s-channel eigenvalue becomes zero at = 0.135, 



and the (ij.2_y2_|_2ia;y-channel eigenvalue becomes zero at 
Xd±2 = 0.35. The s-channel instability corresponds to 
the Fermi surface collapse perpendicular to the dipole 
orientation, and the d-channel one is the biaxial nematic 
instability of the Fermi surface as studied in Ref."^" (If we 
truncate the matrix i^f/Pg at Z = 4, \d±2 '^i^l become 0.95, 
which is essentially the finding in Ref.'^". Therefore, the 
d-channel instability actually occurs earlier than what 
they expect in Refi^) The s-channel instability occurs 
before the d-channel instability with the purely dipolar 
interaction. Nevertheless, the s-channel instability can 
be cured by introducing a positive non-dipolar short- 
range s-wave scattering potential Voo;Oj which adds to 
the Landau parameter of Fq^^q without affecting other 
channels. 

Compared to the compressibility calculation by the 
variational method in Refi^, in which it shows a stability 
condition equivalent to A w 0.42. This is larger than our 
(86) criterion A^ = 0.135 mentioned above. However, the cal- 
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culation in Ref.-- did not take into account the Hartree 
contribution to the ground state energy. Ahhough it is 
zero for rigorously homogeneous systems, it is actuahy 
singular due to the singularity of the Fourier transform 
of the 3D dipolar interaction Vsoiq) as <f — ^ 0. Assum- 
ing an infinitesimal density wave vector q in the ground 
state, the Hartree self-energy contribution appears in 
Foo:o{q)^^ , and becomes most negative at 9q = ^. In 
other words, Refi^ sets F^Q^Q^q) — which overestimates 
the stability of the 3D dipolar gas. Our result is sup- 
ported by the numerical calculation given in Ref*^, in 
which the Fermi surface instability manifests as the on- 
set of an unstable collective mode at As ~ 0.14. For the 
multi-component case, the critical value should be fur- 
ther suppressed by a factor of 1/N because DOS is N 
times large, where N is the number of components. 

Similar static instability also occurs in the 2D dipolar 
system. Let us consider the specific case 9q = Oq, which 
corresponds to the most anisotropic dipolar gas in 2D 
(Eq. [5]) . We numerically diagonalize the Kf-i^^, matrix in 
Eq. [Ml The result shows that an instability would occur 
at A ~ 0.15. 



VI. COLLECTIVE EXCITATIONS IN THE 
DENSITY AND SPIN CHANNELS IN 3D 

In this section, we study the density-density response 
of the dipolar Fermi liquid and the corresponding collec- 
tive excitations. The spin channel collective modes will 
also be studied in multi-component systems. Due to the 
anisotropic nature of the dipolar interaction, the response 
function exhibits anisotropic features. It shows the col- 
lective excitation of the zero sound which only propagates 
within a certain range of directions with anisotropic dis- 
persion relations but become damped in other directions. 

In the following, we will first present the generalized 
dynamical response of the dipolar system for 2 and 3D in 
Sect. IVI Al which is followed by the 3D collective excita- 
tions in Sect. IVIBI In order to obtain a clearer picture 
about the contribution of each mode towards the zero 
sound, we first consider the simplest s-wave channel in 
3D only, where we will see how the zero sound propaga- 
tion is restricted in certain propagation direction relative 
to the external electric field orientation. We then pro- 
ceed to a more quantitative calculation by considering 
the correction from the p-channel. After that, we discuss 
the possible spin collective spin mode in Sect. IVI CI 



same except 61/20 = 6n2Q — ^^^-20 in the channel of {I, m) = 
(2, 0). refers to the equilibrium Fermi surface nematic 
distortion. 

To start with, we consider the standard Boltzmann 
transport equation for the Fermi liquid^''^'' 



dn,- 



UJ - vi^ ■ q deHF,p 



where vi^ is the Fermi velocity in Eq. [25l eHF,p is the 
Hartree- Fock single particle spectrum and f^^ip,]?) be- 
ing the 3D Landau interaction in Eq. 1351 This equation 
is equivalent to 



5v{Q.p) — 



-fpp5v{np,) = 0. (89) 



Lo — vp ■ q J 47r 

After the spherical harmonic decomposition, we arrive 
at the generalized transport equation: 



E 

I'm' 

where 



^^U'Smni' 



E 



xl9,,ram'{^.<l)Ff,Pv.,r.'yyi',n' (90) 



^{2l + l){2l' + 1) 



dn, 



iv^^(%) 



vl'^.q 



U! - 



iQ^.q 



Yl^rrA^p), (91) 



where Fjjj^^ is the Landau parameters defined in Eq. 
[39l Due to the dipolar anisotropy and the propagation 
direction of q, 5vim in different channels are coupled. The 
dispersion of the collective modes can be obtained by 
equating the determinant of the above matrix equation 
to zero. The formalism above is not restricted to dipolar 
system and can by applied to any 3D Fermi liquid with 
azimuthal but not rotational symmetry. 

Generally speaking, in order to calculate the zero 
sound propagation of these anisotropic systems, one has 
to deal with the above infinite matrix equation. However, 
in this 3D dipolar system, it turns out that the physics of 
zero sound is well captured by considering only the s and 
longitudinal p channels. In order to avoid complication, 
in the following, we start by examining the s-channel 
only, where the anisotropic feature of zero sound already 
appears in terms of a limited propagation direction and 
an anisotropic sound velocity. Afterwards, we consider 
the effect of the longitudinal p-wave mode, where a quan- 
titative zero sound velocity is obtained as a function of 
propagation angle and this is in good agreement with a 
numerical study performed in Refi^. 



A. Generalized dynamical response functions 

For the purpose to study collective modes, we define 
Svp as the variation of Fermi distribution respect to the 
equilibrium Fermi surface with the dipolar interaction 
strength A. Compared to the definition Suj; in Eq. 1331 
and (5n(57fc) in Eq. 1711 which refer to the non-interacting 
Fermi surface, their spheric harmonic components are the 



B. The 3D density channel collective mode: zero 
sound 

1. The s-wave channel 

Let us warm up by only keeping the s-wave channel 
component in Eq. 1911 To the lowest order of A, we only 
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keep the anisotropic Landau parameter of -Foo:0 (<?) which 
explicitly depends on the direction of q, and completely 
neglect the anisotropic DOS and Fermi velocity. Thus 
the XQ^gQ(u;,g) is simply the standard textbook result 
which reads at small values of s = uj/vj:^ q ^ 1 as 



xl^Mq,co)^Ni^{l-'-ln^ 



1-s 



+ z|se(s< 1)}. 



(92) 



The collective excitation is determined by the pole of it: 

l + ^oo^o(%)Xoo% = 0- (93) 

s > 1 is needed to ensure the collective mode under- 
damped. 

Eq. [93l gives rise to an anisotropic zero sound velocities 
which explicitly depends on the propagation direction of 
9q as 



where 



4d( 



(94) 



2e I'AP^Ue,); 87rAF2(cos6lg) < 1 
'^AP2(cos6',); SttAFs (cos 6I5) > 1 



,,315 

(95) 

for 0q < 9o or 9q > n — 9o, where ^0 is defined in Eq. 
[31 The zero sound is well-defined around the north and 
south poles of the Fermi surface, where the interaction is 
most repulsive. On the other hand, for ^0 < 9q < n — 9q, 
Foo;o{^q) becomes negative, the solution shows s < 1 
which means the dispersion goes into the particle-hone 
continuum and ceases to be a sharp collective mode. 

In order to refine the above calculation, one has to take 
into account the anisotropic Hartree-Fock single particle 
spectra and thus the anisotropic Fermi surface. In this 
sense, the full Xoo'oo('^i function is given by 



dflp N^^inp){q-\7peHF{p)) 



q- ^p^HFip) + 



(96) 



Nq [1 -I- (cos 9p)] defined in Eq. M 



where N^^^Qp 

and €HF{k) is the Hartree-Fock single particle spectra. 
The angular form factor /(r2p,r2g) is defined as 



q- ^pf-HF{p) 



3 

cos 



AP2(cos6'p)] 



, COSf 



Att sin 29p (sin 9q cos 9p cos 4>p 



(97) 
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FIG. 8: Dispersions of the zero sound s(Qq) = uj{6q)/vjj^ q for 
the pure dipolar interaction including the coupling between 
the s-wave and the longitudinal p- wave modes for A = 
The s-wave approximation is included for comparison which 
is vahd in the limit of a large component number A'' with the 
replacement of A — ^ A^A. When the sound speed hits the par- 
ticle hole continuum, the sound becomes damped. These re- 
sults are in good agreement with a numerical study in Ref.— . 
(The parameters are chosen for a better comparison with the 
Ref.~, where their parameters D = Stt^A are chosen to be 1 
and 3.) 



where the propagation direction q is chosen in the xz- 
plane with the polar angle 9q. Eq. [Mjis solved numeri- 
cally and plotted in Fig. |S| along with the edge of particle- 
hole continuum and the further refined result including 
the p-wave channel correction. (We postpone the discus- 
sion of it until we include the p-channel correction in the 
following subsection.) 

So far, we neglect the coupling between the s-wave and 
other high partial channel Fermi surface excitations. This 
approximation is valid if the Landau matrix elements in 
other channels are small compared to that in the s-wave 
channel. This is justified in the multi-component case 
with a large N in which PboiO ha-s a large A^-enhancement 
compared to all of other matrix elements. Thus the re- 
sult of Eq. inni is correct in the large TV-limit with the 
replacement of A by A^A. However, for the single com- 
ponent dipolar interaction, i.e., A^ = 1, the zero sound 
dispersion is modified from the coupling to the longitu- 
dinal p-wave channel modes as explained below. 
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2. The correction from the coupling to the longitudinal 
p-wave channel 

The coupling between the s-wave and other high par- 
tial wave channel Fermi surface excitations, say, the p- 
wave longitudinal excitation, can significantly change the 
zero sound dispersion. This effect is particularly impor- 
tant if the Landau interaction in the Fi -channel is not 
small compared to Fq. For example, the zero sound ve- 
locity is measured as s = 3.6 ± 0.01 in '^He at 0.28 atm, 
while Fq — 10.8 only gives s — 2. The inclusion of the 
coupling to the p-wave channel — 6.3 gives rise to 
the accurate value of s = 3.6^. Below we consider this 
coupling for the single component dipolar Fermi gas. 

Now, we consider the coupling between the s-wave and 
the longitudinal p-wave channel modes in our 3D dipolar 
system. The formalism developed in Subsect. IVI Al have 
3 p-wave modes; I — 1 and m = 0,±1, which mixes the 
longitudinal and transverse modes. It turns out that it 
is better to reformulate the generalized response function 
by putting the spherical harmonic expansion Yim relative 
to the q direction. Due to the explicit anisotropy, the 
longitudinal and transverse p-wave components should 
be mixed. Nevertheless, the transverse p-wave channel 
mode is overdamped for small positive value of Landau 
parameters, thus they do not affect the zero sound much. 
We only keep the mixing between s-wave and the longi- 
tudinal p-wave modes. 

The corresponding transport equation becomes a 2 x 2 
matrix equation and the collective mode can be solved 
from: 



Bet[l + N{q,u)] ^0, 
where the matrix kernel of N(q,uj) reads 

N(qi.)^( ^o^oo('7',^)«(9l X?o%(?,^)^ 



(99) 



Fi(^Q (q) is the longitudinal p-wave Landau parameter de- 
fined as 



cos^ ^.i^i'i^^^o + sin' O.F^i'^m^^i 



IStt 



= i^lt™=0^2(cOS^,) = _P2(COS0,).(1OO) 

The relevant response functions are correspondingly 
modified as: 



4tt s-f{np,ng) 
47r N§D s-f{np,ng) 



where g*lies in the xz-plane with the polar angle 9q, and 
q ■ k = sin 9q sin 9k cos 4>k + cos 9q cos 9k . 

Before we solve Eq. [98l we can get a qualitative picture 
by approximating the Fermi surface and density of state 
to be spherical first. Under this approximation, in the 
long wavelength limit if — > 0, Eq. [35] can be analytically 
simplified to: 



1+- 



^ln|l±i| 
2 ' 1 - s ' 



(102) 



(101) 



which resembles the usual zero sound condition by con- 
sidering only Fq and Fi parameters in an ordinary 
isotropic Fermi liquid'^^. The difference is that the Lan- 
dau parameters now have an explicit q dependence. Be- 
cause FQo-o(q) and Fn-^oiq) have the same angular depen- 
dence, including the longitudinal p-wave channel mode 
does not change the propagation angular regime, but en- 
hance the zero sound velocity. 

Now, we solve Eq. [98l numerically by taking into ac- 
count the anisotropic fermion spectra and Fermi surface. 
The zero sound propagation as a function of 9q is plot- 
ted in Fig. m The calculation that only involve the 
s-channel is also included for comparison. There are two 
features: a restricted zero sound propagation angle and 
an anisotropic propagation velocity, s. The first feature 
is due to the fact that the sound enters the particle-hole 
continuum for large 9q and is thus damped. It is clear 
that, in terms of the critical angle where zero sound ter- 
minate, the s-wave approximation is qualitatively justi- 
fied. However, when we consider the zero sound speed as 
a function of 9q, the longitudinal p-wave mode modifies it 
considerably. These results are in a good agreement with 
a fully numerical calculation based on the same Boltz- 
mann transport theory'^^ . This indicates that the physics 
of collective excitation of the 3D dipolar system is well 
captured by considering only the s and longitudinal p- 
channels. 

The zero sound velocity is affected by both the value 
of Landau parameters and the Fermi velocity as shown 
in Eq. For the relatively large value of A = the 
Landau parameters play a more important role. Thus 
the sound velocity is largest for q along the north and 
south poles where -Fo(k)('L) ^-^d Fivoi^ largest. As 
the propagation deviates from them, the sound veloc- 
ity first becomes softened. As the velocity hits the up- 
per edge of the particle-hole continuum, the zero sound 
ceases to propagate. At a smaU value of A = Fig. [5] 
shows the upturn of zero sound velocity before it hits the 
particle-hole continuum, which is largely determined by 
the anisotropic Fermi velocity. We also plot the eigen- 
vector for the zero sound mode for the solution of Eq. 
1981 The s and the longitudinal p-wave components are 
comparable to each other. 
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FIG. 9: The ratio between the longitudinal p-wave and the 
s components in the eigenvectors of the 3D zero sound, i.e. 
R = Svio-o/Svoofi, for A = and ^i^-. 



FIG. 10; 3D collective spin mode in the longitudinal p-channel 
for A = The sound speed continuously decreases as 9q 
increases and terminates at certain critical angle. 



C. 3D Collective spin mode in the longitudinal 
p-wave channel 

In addition to the zero sound mode, collective exci- 
tations may also exist in the spin channel. Because 
-^(KH)" ~ ^ (^^^ definition in Eq. [50)) . there is no well- 
defined s-wave channel spin excitations. Nevertheless, 
the longitudinal p-wave channel mode becomes a propa- 
gating mode. The formalism to determine its dispersion 
is very similar to Eq. [90] by replacing F^,.^ to F^,.^. 
Only keeping the longitudinal p-wave channel, we have 

l + ^n1'd"(9)x?f?oo(9,^) = 0. (103) 

^ii^o'^i'l) — -^n;O^2(cos0^) which is the same as F^(^Q^{q) 
at the Hartree-Fock level, and Xi?oo ^^^^ same as 
in Eq. IIOII The dispersion for this p-wave longitudinal 
spin collective mode is plotted in Fig. [10] for A = 
Again, we have an anisotropic sound velocity subject to a 
finite propagation regime, where the excitation is damped 
when its dispersion enters the particle-hole continuum. 

In A'^-component dipolar gases with the SU{N) sym- 
metry, there are iV^ — 1 branches of longitudinal p-wave 
spin model. They are essentially the fluctuating longitu- 
dinal spin current mode. The is a novel collective mode 
which has not been observed in condensed matter sys- 
tems before. 



VII. CONCLUSIONS AND OUTLOOK 

In summary, we studied the anisotropic Fermi liquid 
states of the cold atomic dipolar Fermi gases which pos- 
sesses many exotic features. The feature of the ciz2_3,.2 
symmetry of the 3D dipolar interaction and the d^2_y2 
symmetry of the 2D dipolar interaction nicely render 
most analytic calculations possible. The same anisotropy 



exhibits at the leading order of perturbation theory in 
many ways such as the Fermi surface anisotropy, Fermi 
velocity and effective masses. The Landau parameters in 
isotropic Fermi liquid states are generalized to the Lan- 
dau interaction matrix which has the tri-diagonal struc- 
ture as a result of the Wigner-Eckart theory. Physical 
susceptibilities receive renormalization from the Landau 
interaction matrix. With large dipolar interactions, dipo- 
lar Fermi surfaces become collapse along the directions 
perpendicular to the dipolar orientation. We also studied 
the collective mode in the density and spin channels. The 
zero sound exhibits anisotropic dispersion relation with 
largest propagation velocity along polar direction. The 
p-wave longitudinal spin channel mode is a well defined 
propagating mode for any propagation directions. 

So far, we have only considered the spatial homoge- 
neous systems. However, the realistic experimental sys- 
tems have confining traps, which should bring corrections 
to properties discussed above as studied in Refs^^i^i^^. 
In the case of soft confining trap potentials, we can treat 
the inhomogeneity by using the local density approxima- 
tion. Around each position r , we can define a local Fermi 
energy and the local Fermi surface which are determined 
by the local molecule density. The molecule density is the 
highest at the center of the trap, thus both the absolute 
interaction energy scale and the dimensionless interaction 
strength are strongest at the center. Thus the thermody- 
namic instabilities are strongest at the center of the trap 
and become weaker as moving towards the edge. For 
the collective excitation of zero sound, the local sound 
velocity is largest in the center. As a result, the local 
propagation wavevector increases from the center to the 
edge to maintain the excitation eigen-frequency the same 
in the entire trap. The collective modes have more server 
damping in the edge area because the molecule density is 
low and thus is less quantum degenerate than the center. 
A detailed study of all the above effects will be presented 
in a later publication. 
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Other open questions to be studied in the future in- 
clude the issues of higher-order (i.e. beyond Hartree- 
Fock theory) interaction corrections and finite tempera- 
ture corrections. The theory developed in this work is 
strictly valid for weak interaction and low temperatures. 
For example, an important aspect of the dipolar Fermi 
gas is the finite life time and the Fermi liquid wavefunc- 
tion renormalization factor Z of quasiparticles. Both of 
them naturally also develop anisotropies due to the dipo- 
lar interaction. However, they begin to appear at the 
level of the second order perturbation theory beyond the 
Hartree-Fock level, thus their anisotropies should be of 
an even higher order spheric harmonics than the d-wave. 
Another issue is the screening effect in the dipolar Fermi 
gases which may change its long range negative in 3D. 



These effects will be deferred to another paper—. 
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